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In this paper we solve two equivalent time optimal control problems. On one hand, we design 
the control field to implement in minimum time the SWAP (or equivalent) operator on a two-level 
system, assuming that it interacts with an additional, uncontrollable, two-level system. On the 
other hand, we synthesize the SWAP operator simultaneously, in minimum time, on a pair of two- 
level systems subject to opposite drifts. We assume that it is possible to perform three independent 
control actions, and that the total control strength is bounded. These controls either affect the 
dynamics of the target system, under the first perspective, or, simultaneously, the dynamics of both 
systems, in the second view. We obtain our results by using techniques of geometric control theory 
on Lie groups. In particular, we apply the Pontryagin Maximum Principle, and provide a complete 
characterization of singular and non-singular extremals. Our analysis shows that the problem can 
be formulated as the motion of a material point in a central force, a well known system in classical 
mechanics. Although we focus on obtaining the SWAP operator, many of the ideas and techniques 
developed in this work apply to the time optimal implementation of an arbitrary unitary operator. 
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I. INTRODUCTION 


A. Statement of the problem 


In this work we consider the controlled dynamics of a pair of interacting two-level quantum 
systems, where the purpose of the control action is the generation in minimum time of a 
desired unitary operation on one of these systems, which in the following will be called the 
target system (T). The other system represents an unavoidable disturbance or environment 
(E). The dynamics of the whole system is given by the Schrodinger operator equation 

X{t) =-iH{u)X{t), A(0) = /®/, (1) 

where X{t) describes the evolution of the density matrix of the composite system from its 
initial to its final configuration, p(0) p{t) = X{t)p{0)X\t), and I is the identity operator 
for two-level systems. The initial state is assumed to be factorized, p(0) = Pe{0) 0 Pt(0), 
and H{u) is the Hamiltonian of the system, given by 


H{u) = / 0 {uxSx + UySy UzSz) + 2ojqSz 0 Sz. 


( 2 ) 


The first term in H{u) contains three possibly time varying control actions Uk = Uk{t), 
k = x,y, z, with L 2 norm bounded in strength, ul + Uy + uj, ^ 7 ^. These control actions are 
associated with the generators Sk = \(Jk, k = x,y,z, where ak are the Pauli matrices 


^ T. 



ay 


0 -i 
i 0 


(Tz = 


1 0 
0 -1 


( 3 ) 


with commuting relations [ak, ai] = 2iam, and {k, I, m) is a cyclic permutation of (x, y, 2 :) j29|. 
The second term in (|2]) represents an Ising interaction between the two systems; we assume 
Uo 7 ^ 0 , otherwise the problem reduces to that of controllii^ an isolated two-level system with 
bounded control, which has been solved in previous work [l|,l2|. Without loss of generality we 
can take Uq > 0, as we shall shortly see. We neglect the free evolution of the environmental 
two-level system. It is meant that terms on the right of the tensor product are associated to 
the target system. Because of the block structure of the Hamiltonian (j2]) and of the initial 
state A( 0 ), the evolution operator has a block structure at any time t, 


X(t) = diag(Xi(t),X 2 (t)) = Pi0Xi(f) + P 2 0 X 2 (t), (4) 
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where the projectors Pi, P 2 are given by 


= + P2 = ^(/-(t^), (5) 

and the dynamics of the two blocks is given by 

Xi(t) = —iioj^Sz-\-UxSxUySyUzSz)Xi(t), Xi(0) =/, 

X2{t) = -i{-uoS^ + UxSx^-UySy + UzSz)X2{t), X2(0) = J. (6) 

The dynamics of the target system is obtained by taking the partial trace of p{t) over the 
degrees of freedom of the environmental system; since TrPjPj = 6ij, we hnd 

PT{t) = TTE{p{t)) = TTE{PlPE{0))X,{t)pT{0)Xl{t) + Tte {P2Pe{ 0)) X2{t) pT{0)Xl{t). ( 7 ) 

Our task is to hnd the minimum time tf such that this evolution reproduces a given unitary 
operator Xf E SU{2) acting on the target system, 

PT{tf) = XfPT{Q)Xy ( 8 ) 

In principle, the initial state of the environmental qubit p_b( 0) is unknown. Therefore, Eq.s 
d?]) and dH]) must hold for every Pe( 0). A necessary and sufficient condition for this is 
that Xi{tf) = and X 2 {tf) = ^Xf, with all possible combinations of signs. In view 
of the evolution equations dS]), we conclude that the problem of optimally steering a two- 
level system in interaction with a second two level system is equivalent to the problem of 
optimally driving two independent two-level systems, subject to opposite drifts, through the 
same control actions. Moreover, we see that our choice Wq > 0 is not restrictive, since the 
case Wo < 0 amounts to interchanging the two systems 1 and 2. 


B. Motivation 


The application to the study of quantum mechanical systems of methods of geometric 
control theory, and in particular optimal control theory, is a longstanding practice. These 
methods are particularly suitable in applications to the synthesis of specific operations in 
quantum information processing j^, the control of nuclear spins in nuclear magnetic reso¬ 


nance 


or the optimization of specific reactions in atomic or molecular physics. 
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In this context, two-level systems are often the bnilding blocks of larger strnctnres. There¬ 
fore, it is of great relevance to understand to what extent it is possible to manipulate them, 
especially when they are immersed in a dissipative environment (see 5l4l3| and references 
therein for some applications of optimal control theory to quantum systems). The standard 
approach for their control is given by dynamical decoupling techniques [l^, which, however, 
usually require unbounded control actions. Whenever a bound on the controls cannot be 
neglected, optimal control represents an alternative tool to dynamical decoupling. 

A dissipative environment is usually modeled by a thermal bath with inhnitely many 
degrees of freedom. However, for some specihc systems, it is meaningful to model the 
external world by means of a hnite-dimensional system. A noticeable example is provided 
by the study of NV centers in diamonds, which are currently of great interest for applications 
in quantum information technology (see, e.g. [^, 16|). The control task is to drive a central 
spin independently of the environmental spins, whose coupling strengths to the central spin 
generally depend on their distances from it. The problem considered in this work is the 
simplest instance of this situation where only one environmental spin exists. See [l^ for 
a numerical analysis of the dependence of the optimal time of transition on the number of 
environmental spins. 

In all these scenarios it is desirable that a given state transfer occurs in minimum time. 
This is both to increase the efficiency of quantum operations which typically require a cascade 
of elementary gates and to minimize the effect of an (un-modeled) environment. 

Besides applications in quantum mechanics, the optimal control of SU{2) operations is 
strictly related to the optimal control on SO{3), being SU{2) the double cover of SO{3). 
Since SO{3) is the group modeling the attitude of a rigid body, our analysis is also of general 
interest in the study of classical mechanical systems. 


C. Plan of the paper 

We are able to derive analytical expressions for the time optimal control strategies, and 
the corresponding optimal times, for target operations equivalent to the SWAP operation, 
in the terms described in Proposition III. 11 and for a specihc range of values of uo and 
7 . Nonetheless, whenever appropriate, we present partial result valid for general target 
operations Xf. 
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In Section mi we discuss the Lie algebraic structure of the probleuu and briefly review 


the Pontryagin Maximum Principle (PMP) of optimal control theory 1^, ll9|] in the context 


of control on Lie groups. This fundamental result provides the theoretical framework for 
our investigation. The controls satisfying the PMP are candidate optimal and they are 
called extremals, and they are classihed in nonsingular and singular extremals. Usually, the 
procedure for solving an optimal control problem by means of the PMP consists of several 
steps, which in our specihc case are given by: (i) solution of the dynamical equations derived 
from the PMP; (ii) determination of the extremals; (hi) solution of the system dynamics 
for these control strategies, and computation of the corresponding time of transition to 
the desired target operator; (iv) determination of the optimal strategy by comparison of 
the different extremal strategies. In general, some of these steps cannot be carried out 
analytically, and one has to resort to numerical computations. In the present case, numerical 
analysis and the structure of the problem will suggest a conjecture, which, when adopted, 
allows an analytical result of the optimal control problem when Xf = Xswap- 

We shall present a complete investigation of the aforementioned points (i) - (iv) for 
the system of interest here. This process is complicated by the presence of singular arcs. 


which require a separate analysis (see 20 l-l 22 [| for applications of optimal control theory 


to quantum mechanical systems where singular solutions exist). In Section IIIII we show 
that, in both cases of nonsingular and singular extremals, the basic equations derived by 
applying the PMP can be interpreted as the motion of a hctitious material point subject 
to a velocity-dependent force. If we choose the hnal target operator as Xf = Xswap, this 
force becomes conservative, and the system is equivalent to an integrable classical central 
force problem. This fact provides a parallel between this optimal control problem and the 
(deeply investigated) analysis of the motion of a material point in a central force 23 1 . 

At this point, we hnd convenient to limit our attention to the case Xf = Xswap, and 
show that the extremal trajectories correspond to motions with null angular momentum in 
the effective central-force problem. This result is valid for both cases of nonsingular and 
singular extremals, which are separately investigated in Sections IIVI and |Vl A complete 
treatment of singular extremals can be performed for arbitrary values of 7 and Uq, but 
the investigation of nonsingular extremals is increasingly difficult as ^ decreases, therefore 
we limit our attention to cases where analytical solutions are possible. In Section IVIl by 
comparing the different types of extremals previously found, we solve the optimal control 
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problem for X/ = Xswap in the aforementioned range of values of wq and 7 , and provide 
the optimal control strategies and the optimal time. In Section fyill we discuss our hndings, 
relate them to existing results, and conclude. 

For sake of clarity, the most technical results are reported in Appendices. Also, the most 
signihcant results, which are frequently referred to throughout the paper, are presented as 
Propositions, and formally proved. The reader not interested in the computational details 
could directly move to Section IVTl where the main results are summarized. The investigation 
of cases other than X/ = Xswap will be presented in a forthcoming paper, complementing 
the analysis presented here. 


II. THE LIE ALGEBRAIC STRUCTURE OF THE PROBLEM 


As anticipated in the Introduction, the optimal control problem that we consider in this 
work is the following: given the dynamical system ([ 6 ]) and a target operation X/ G SU{2), 
we want to hnd a control strategy Uk(t), k = x,y, z, such that Xi(tj) = ±Xj and X 2 (tf) = 
±Xi{tf), and tf is minimum 


30| . First of all, we prove that, for every 0 G M, the operator 


Xj = 


( 9 ) 


is reached in the same optimal time as Xf. 

Proposition II. 1. Consider a pair of operators Xf and Xf G SU{2) which satisfy ([2]) for 
some (p. Assume that there is a control strategy Uk = Uk{t) {k = x, y, z), with u^ + Uy + uj. ^ 
7 ^, such that Xi and X 2 are mapped to Xf (or its negative) in time tf. Then there is a 
control strategy Uk = Uk(t) {k = x, y, z), with ul + Uy + ul ^ 7 ^, such that Xi and X 2 are 
mapped to Xj (or its negative) in the same time tf. 

Proof: If we dehne Xi = and X 2 = we can write 

Xi(f) = -i{uJoSz + UxSx + UySy + UzSz)Xi{t), Xi(0) = /, 

^2{t) = —'i{—^oSz + UxSx-\-UySypUzSz)X2{f), X 2 ( 0 ) = /, ( 10 ) 

where Ux = UxCos^p + Uy sin (f, Uy = —M^-sincp + Uy cos (p, and Uz = Uz. These controls 
satisfy the required constraint since u^. + = ul. + + ul, therefore they represent 

an admissible control strategy. Moreover, the dynamics ffTOj) is formally the same as (jO]), 
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therefore the control strategy Uk = Ukit) generates the desired transformation Xi{tf) = 
and X 2 {tf) = ±Xi{tf). 


This result is a consequence of a symmetry of the problem, dehned by the transformations 
dfc = for k = x,y, and = a^- This change of representation of the Pauli 

matrices does not affect the system, since it merely corresponds to a redefinition of the 
controls, but it leaves invariant the equations of motion ([6]). The particular form of the 
controls in the dynamics is crucial for the existence of this symmetry, which is preserved 
even if there are only two independent controls affecting and ay (that is, = 0). 

Following Proposition III. 11 we will consider equivalent two target operators which satisfy 
(P|) for some (p. Therefore, it is sufficient to solve the optimal control problem for one of 
them to obtain a complete description of the optimal control strategies and optimal time for 
all of them. We will rely on this fact, when we will fully solve the case with Xf = Xswap 
(or equivalent operators). In the representation induced by ([3]), operators related by ([9]) 
differ by a phase in the off diagonal terms. 


It is well known in quantum control theory 2^ that information on the structure of a 


control problem can be obtained by analyzing the Lie algebra associated to the system. We 
dehne the block-diagonal matrices: 


Bh — 



C,= 



( 11 ) 


where Sk = ^ak, and 0 is the 2x2 null matrix, and k = x,y, z. In terms of them, the 
evolution ([1]) takes the form 


X{t) = -i{uoC^ + + UyBy U:,B^)X{t), X{Q) = I ® I. 


( 12 ) 


Since we want to hnd the controls steering in minimum time tf the evolution operator to 
X{tf) = dia.g{Xf, Xf), or X{tf) = diag(Xj, —Xj), or their negatives, where Xf G SU{2), 
we hrst apply results on the controllability of the system, to answer the question of whether 
this type of state transfer is possible. This is indeed the case since the Lie algebra generated 
by the Hamiltonians Cz, B^, By and Bz in (IT^ is, in the considered representation, the Lie 
algebra of block diagonal 4x4 matrices with arbitrary 2x2 blocks in su(2). The set of 
reachable operators is the associated Lie group of block diagonal 4x4 matrices with 2x2 
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blocks in SU{2). This group contains the desired hnal conditions for any special unitary 
operator Xf. Therefore, system flT^ is controllable on the Lie group of operators of the 
form diag(Y,Z), where Y, Z E SU{2). This group is compact, semisimple, and isomorphic 
to SU{2) © SU{2). The associated Lie algebra, isomorphic to su(2) ©su(2), is given by 
1 = b © c, where 

b = span{Bj;, By, B^), c = span(C',j,, Cy, C^) (13) 

satisfy the commutation relations of a Cartan decomposition of 1, 

[b,b]Cb, [c,c]Cb, [b,c]Cc. (14) 

More specihcally, the commutation relations of the operators in (ITT]) are given by 

[Bk,Bi] = iBm, [Ck,Ci] = iBm, [Bk,Ci] = iCm, (15) 


and (/c, I, m) is a cyclic permutation of (x, y, z). The controllability properties of the system 
along with the fact that the set of possible values for the control is compact allows us to use 
the standard Filippov’s existence result for time optimal control problems (see, e.g., Hi to 
conclude that the time optimal control exists for every hnal condition in SU (2) © SU (2) 31 |. 

To solve our problem, we shall apply the classical neces sary conditions of optimality given 
by the Pontryagin Maximum Principle (see, e.g., and 26, 271 for application to quantum 
systems on Lie groups). In order to state them in a form appropriate for our goals, we set 
up some dehnitions. Given M G 1 (usually called the costate), and X a solution of flT^ . we 
dehne the time-dependent coefficients 

bk = i{M, X^BkX ), Cfc = ^(M, X^CkX ), (16) 

where (^4, B) = Tr(ili3l), and k = x,y, z. 


Definition 1. The Pontryagin Hamiltonian is dehned as 

H{M,X,v^,Vy,V:,) = uqC:,+ ^ Vkbk. (17) 

k=x,y,z 

Theorem 1. - Pontryagin Maximum Principle (PMP) - Assume that the evolution of X is 
described by [Wt) . and that Ux,Uy,Uz is the time optimal eontrol strategy steering X from 
the identity I ® I to diag(X/,X/) or diag{Xf,—Xf). Denote by X{t) the eorresponding 
trajectory. Then there exists an operator M E I, M ^ 0, a constant X, and functions bk, 
Ck as in m. such that H{M,X,Ux,Uy,Uz) ^ H{M,X,Vx,Vy,Vz) for every Vx, Vy, Vz such 
that vl + Vy + vl ^ 7 ^. Moreover, the Pontryagin Hamiltonian is constant and egual to X. 







Controls satisfying the conditions of Theorem [T] are called extremals, and they are can¬ 
didate optimal controls, since the PMP is a necessary condition for optimality. They are 
called normal extremals if A 7 ^ 0, abnormal otherwise. The arcs of the extremal trajectories 
where the Pontryagin Hamiltonian does not explicitly depend on the controls are called 
singular, and in this specihc case they are described hy = by = = Q. In general an 

extremal trajectory can contain both singular and nonsingular arcs. We shall call singular 
a trajectory which contains at least one singular arc, nonsingular otherwise. 

The procedure to hnd the time optimal control consists of computing all the extremals, 
and then compare the values of the time needed to reach the desired hnal condition and 
choose the control which gives the minimum time. 

For both singular and nonsingular arcs, and, in fact, for every control, the functions 
dehned in ffT 6 |) satisfy the following differential equations, 

bx "^zby T '^ybz-i 

by UqCx T u^bx Uxbz, 

bz 'Uxby Uybx, 

Cx ^oby T UyCz "^z^yi 

Cy ^obx '^xCz T '^zCx: 

Cz ’^x^y ^y^X' 


They can be derived by differentiating ffTbl) . taking into account the evolution flT^ . and the 
commuting relations flTSl) . 

For some given initial conditions 6^(0), and Cfc(O), k = x,y,z the hnal conditions bk{tf) 
and Ckitf), are obtained by taking Xf = X{tf) in flT^ . They read 


bk{tf) — ''^^Xkibi{0), Ck{tf) — (19) 

i i 

where X G SO{3) is a matrix representation of the element associated with Xf. This follows 
from the fact that SU{2) is the double covering group of 50(3). 

It turns out that fll7p is not the only constant of motion on extremal trajectories. 


Proposition II.2. System flTH]) admits the following integrals of motion 32l |: 


^=Y (bl + cl), L = uo ^kCk- 

k=x,y,z k=x,y,z 


( 20 ) 
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Proof: The fact that E = L = 0 is a. direct consequence of fflSj) . 


In the next section we will provide a useful interpretation of these constants. To complete 
this section, we derive some results which will be needed in the following. First of all, we 
characterize the structure of extremal control strategies. On nonsingular arcs, since the 
Pontryagin Hamiltonian is linear in v^, Vy and Vz, its maximum in the set v^-l-Vy + vl ^ 7 ^ 
must be on the border of this set. By applying the Lagrange multipliers method we hnd 
that the extremal controls satisfy 


Uk = 7 - 




k = x,y,z 


( 21 ) 


where jj^b = + by + b^. On nonsingular arcs, fib can vanish only in isolated points. From 

standard theorems on ordinary differential equations, it follows that bk must be continuous, 
and then the controls Uk must be piecewise continuous. 


Remark II. 3 . At the beginning of this section, we have proved that two target operators Xf 
and Xf satisfying (I9]) are reached in the same minimum time. One possible interpretation of 
this result is to think of (|9]) as a symmetry transformation which preserves the Lie algebraic 
structure of the problem. Accordingly, if the matrix M G 1 appearing in the PMP is 
associated with Xf, the corresponding operator M G 1 associated with Xf, must be given 
by 


M = (^eM. 


( 22 ) 


The coefficients flTHll are transformed as = bx cos '•p + by sin (p, by = —bx sin p + by cos p and 
bz = bz, where bk, k = x,y, z, are associated with M (the same transformation describes the 
map of the Ck coefficients into Ck), and the Pontryagin Hamiltonian is unchanged. 


III. THE COSTATE DYNAMICS AS A CENTRAL-FORCE PROBLEM 

In this section we reinterpret the costate dynamics, i.e. Eq.s (fT 8 |) . as a central-force 
problem. This parallel provides a useful bridge between our problem and a well know 
scenario in classical mechanics, which has been deeply investigated in the past. While we 
will find this connection especially useful in the special case Xf = Xswap or equivalent 
operator, we believe that it is signihcant for more general hnal conditions. Therefore, for 


10 



sake of completeness, initially we consider a generic target operator, and then we specialize 
our analysis to the case-study of X/ = Xswap- 

Assume we are on a nonsingular extremals. By using the general form of the controls 
(ED, we see from fflS]) that = 0, and then bz = bz{0) is a constant. This implies that 

fib (bxbx T ^y^y') (^ 2 ) 

which is equivalent to the constancy of the Pontryagin Hamiltonian on nonsingular ex¬ 
tremals: 

H{M,X,u^,Uy,u^) = A = a;oC^ + 7/^6 = ^oC^(O) + 7/^b(0). (24) 

This result is embodied in the last equation in flTSl) . The other equations reduce to 

bx UjQCy^ 

by 

7 

Cx (^oby T {byCz bzifi)Cy^ ^ 

7 

Cy ^obx (bx^^z bziO'jCx^ j ( 2 ^) 

hb 

which, on account of fl 2 T|) . can be written as 

bx , = -(wq + + 7^—- 

fJ^b f^b 

by = —(wq + + yA— 762 ( 0 )— (26) 

fib f^b 

These non-linear equations can be interpreted as modeling the motion of a material point 
of unit mass, described by a “position” vector b = {bx,by) = {bcos9,bsm9), in Cartesian 
or polar coordinates, respectively. This point is driven by a central force, the hrst and 
second terms in the right hand side of (l26|) . plus a term dependent on “velocity” (the last 
contribution in both equations). By using a compact notation, 

b = F(b) + f(b,b), (27) 


where 


F(b) = -(wq + 7^)b -h 7A 




f(b,b) ^ -ixh(())ay 




(28) 


The velocity-dependent term is not a dissipative contribution. In fact, it is possible to 
interpret the integral of motion E of the system, given in Proposition III.21 as the associated 
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energy. This can be seen by introdncing the potential associated to the radial part of the 
force, 

m = + 7 ") {b^ + bl(0)) - + b'HO) + iA^ (29) 

with F(b) = —gradi,ff(fo). It follows that 

B = (7(6) + i(i-^ + (6(j)^), (30) 

which jnstihes the interpretation of E as the energy (potential pins kinetic energy). The 
energy is conserved, and there is no dissipation. 

Since we have fonnd bz = bz{0), and by virtne of (IT^ . we can conclnde that b(tf) = 6(0). 
Since E is an integral of motion, we mnst also have c^(t/) = 0 ^( 0 ) (alternatively, this is a 
consequence of the invariance of the Pontryagin Hamiltonian). 


A. The case Xj = Xgwry^^p 


If we assume Xf = Xswap or any equivalent operator parameterized by g) according to 
Proposition III. 11 we have some simplihcations. Therefore, here and in the following sections 
we will limit our attention to this special case, and then comment about the general case 
again in Section IVIIl Since we consider special unitary target operators, we conventionally 
write Xswap = icTy, and then the aforementioned family of equivalent operators is given by 


x,= 




(31) 


where ip is an arbitrary real number, which can be taken in the interval [ 0 , 2n) without loss 
of generality. The SO{3) matrix introduced in 0191) is given by 


^ — cos 2(p sin 2 (y 9 0 ^ 


A' = 


V 


(32) 


sin 2ip cos 2ip 0 
0 0-1 

which is a rotation of angle tt about the axis dehned by the unit vector (sin 93 , coscp, 0). For 
the hnal coordinates we must have 


/ 


bx{tf) = — 6 a;( 0 ) cos 2 (p + 6 j^( 0 ) sin 2 (p, 
by{tf) = 63 ,( 0 ) sin 2 (p + 6 j^( 0 ) cos 2 (p, 

bz{tf) = -6,(0), (33) 
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or, in polar coordinates, b{tf) = 6(0) (as we already know), and 0(tf) = 7r — 0(0) — 2cp. Similar 
relations hold for Ck{tf). Consistency with the former constraints requires bz{tf) = 6^(0) = 0 
and Cz{tf) = Cz{0) = 0. This implies that the Pontryagin Hamiltonian equals A = 76 ( 0 ), the 
velocity-dependent term disappears, f(b,b) = 0, and the central force simplihes to 

F(b) =-(^ 0 + 7^)b + 7Ap (34) 

The associated potential is given by 

U{b) = ^iul + -f^)b'^--fXb, (35) 

which belongs to the class of integrable potentials, with solutions given in terms of elliptic 
integrals. Initial and hnal positions and velocities are related by an orthogonal transforma¬ 
tion, the upper 2x2 diagonal block in (1321) . Moreover, the integral of motion L, given in 
Proposition III.21 takes the form 

L = bH, (36) 

and then it can be interpreted as the angular momentum of the system. 

Because of its particular simplicity, which highly limits the need of numerical analysis, in 
the reminder of this paper we will fully explore this case. As anticipated in the Introduction, 
we now separately consider nonsingular and singular extremals. 

IV. NONSINGULAR EXTREMALS FOR Xf = XswAP 

In this section we limit our attention to nonsingular extremals, for which /if, can vanish 
only at isolated points. Following the analysis of the previous section, we know that it is 
possible to express the costate dynamics in terms of elliptic integrals. Unfortunately, this 
is not enough to provide useful expressions for the integrals of ([ 6 ]), that is, suitable expres¬ 
sions for the comparison of the extremal trajectories. Nonetheless, numerical computations 
suggests that the minimum time trajectories are characterized by L = 0. We are not able 
to prove this result analytically, but we have found numerical evidence of its validity. We 
refer to Appendix for more details on the numerical analysis supporting this fact. 

Conjecture 1. If Xf = Xswap, o,nd an extremal control strategy is optimal, then the 
corresponding trajectory satisfies L = 0. 
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Under this assumption, the dynamics of the costate greatly simplihes, and it is possible 
to compute explicit expressions for the evolution in the Lie group. 

The following argument is based on Proposition III. 31 If we assume that the optimal 
control strategy leading to Xf = Xswap requires b( 0 ) = ( 6 a;( 0 ), &y( 0 )), then an equivalent 
operator, differing from Xf hj a. phase ip in the off-diagonal elements, is associated with the 
counter-clockwise rotation by an angle ip of b(0). In particular, one of these operators has 
to be associated with b(0) = ( 63 ,( 0 ), 0), with 63 ,( 0 ) ^ 0. According to Proposition lII.il there 


is no loss of generality in choosing_ 


equivalent to SWAP will follow 331. 


his operator as Xf: the general analysis of any operator 


Since we have assumed L = 0, we can write b(f) = ( 63 ,( 6 ), 0), and then 6 ^( 6 ) = 0, 
which implies by{tf) = 0. But validity of fl5^ imposes that sin2(/? = 0, which can be 
satished either hy ip = 0 or ip = |. These two cases correspond to X/ = iay and Xf = 
ia^, respectively. Before analyzing them, we determine the structure of extremal control 
strategies for nonsingular arcs. 

Since /Uf, = |6a,|, we can rewrite the hrst equation in fl26l) as 


61 = -(Wj + A)h + sign(6i)7^i>i(0). 


(37) 


with general solution 


63 ,( 6 ) = Aj cos ut + Bj sina ;6 -|- sign( 63 ,) j 63 ,( 0 ), (38) 

where = cuq + 7 ^. The coefficients Aj and Bj, j = 0,1,2,... are determined by the initial 
conditions, and by the requirement that 63,(6) and 63,(6) are continuous functions (whenever 
63, vanishes, new coefficients must be taken into account). The initial pair of coefficients is 
given by 

Tlo = (—)h,(0), B„ = ±i^Vl-6.(0)2, (39) 

\ UO / CO 

and Bq has been written under the assumption that E = ^, since the value of this constant 
does not affect the optimal control problem. The sign of Bq is the same as the sign of 63 ,( 0 ). 
We have two possible solutions fl38l) . depending on the sign of Bq. All the other coefficients 
Aj and Bj are completely determined by the continuity requirements, and their explicit 
expressions are not needed for further developments. Following fl2T|) . the optimal control 
strategy is necessarily bang-bang, with 


M3, = 7 sign (63,), Uy = u^ = 0 . 


(40) 
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If they exist, the switching times are given by the zeros of fl38l) . which depend on 6a;(0). We 
denote these switching times hy ik = i + ki, k = 0,1,2,..., with t ^ i ^ 0. The zeros of 
fl35]) satisfy 


cos ujt± 


-A„C„ ± BoVAJ + - CS 

Al + Bl 


sin ut± 


-BqCq t A,^AI + Bl - Cl 

+ Bl 


(41) 


with Ao and Bq as in fl39|) . and Cq = (^)^^a:(0). We can assume that and t_ are, 
in absolute value, the smaller and next to smaller zeros, not necessarily in this order. If 
&a;(0) 7 ^ 0, by considering both positive or negative values for Bq in fjTTl) . we hnd that 
sina;t+ ^ sina;t_. Moreover, for Bq ^ 0 we hnd cosa;f+ ^ cos ut-, cos ut- ^ 0 and 
sinci;t+ ^ 0, and for Sq < 0, we hnd cosa;f+ ^ cos cutcos out + ^ 0 and sina;t_ ^ 0. From 
these properties we can conclude that, in any case, ^ 0 and ^ 0 , so that 


cosut = cos ut_, smut = smut_ 


(42) 


and t = (t_ — t+), leading to 


cos ut = 


2CS -AS- Bl 

+ Bl 


sin cut = 


2C,^Al + Bl - Cl 
Al + Bl 


(43) 


^0 T -^0 ^0 T -^0 

We don’t need explicit expressions of i and t in terms of wq, 7 and 6a;(0), and Eq.s 
and (ITT|) will be sufficient to fully specify the control strategies in the case of interest 
for this work. In the special case bx{0) = 0, we have Aq = Cq = t) and Bq = ±^, from 
which we can evaluate t = t = -. For further reference, we observe the following: 


Proposition IV. 1. From sin cut ^ 0 we have that - < t < —. Moreover, 0 < t < 


In the following, we separately consider the two aforementioned cases Xf = ioy or Xf = 

A. Case with Xf = iay 

By using (p = 0 in ([33]), we derive the hnal condition bx{tf) = — 63 ,( 0 ). An explicit 
computation proves that, if 63 ,( 0 ) = 0 and the control Ux is constant, it is impossible to 
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reproduce Xf. Therefore, control strategies associated with these nonsingular extremals 
have an odd number of switching times. We adopt the following dehnitions: 

5+ = -( 7 S', + n;o^,), 5_ =-( 7 ^, - (44) 

OJ bJ 

and with uj = -^(Uq + 7 ^. By integrating the dynamics ([ 6 ]), we hnd 

V(t) = Ki(()(y+(riyi(())”y+(«), 

= y 7 f)(y-(f)yt(‘))”y-(<). (45) 

where n = 0 ,l, 2 ,..., and we have used the notation 

y+(t) = y_(t) = (46) 


By requiring that Xi{tf) = —X 2 {tf) = ±Xf, from the pair of equations fITSD . we obtain 


(y+(f)yi(i))" = ±y_(t)wy7(). 


(47) 


and the total time of evolution on these nonsingular trajectories is tj = 2{t+nt). It turns out 
that the other possibility, Xi{tf) = X 2 {tf) = ±W/, is inconsistent. To compactly present 
the treatment of this case, we hnd convenient to dehne two functions a = a{uo,'y,t) and 
(3 = (3 satisfying 


sin a = 


(^o\ . 

— sinwr. 


cos a = 


. 2 - 

— sm wr, 


\ijj J V \UJ J 

1 - 

sin f3 = — sec a cos ur, cos (3 = — sec a sin wr, (48) 

OJ 

where f = ^. These dehnitions are meaningful only when cos a 7 ^ 0, which is always the 
case if 7 7 ^ 0. To prove that they are consistent is a lengthy but standard problem of 
trigonometry. We can now write in compact form 


y+(f)y!(t) = cos2aI + zsin2Q; (cos /3ay + sin /3cr^), 


(49) 


from which it follows that we can express the left-hand side of fl47)) as 

(Y+{t}Yl{F)^ = cos2naI + ism2na{cos/3ay + sin/3az). (50) 

On the other side, the right-hand side of flT7|) is computed by considering (146|) and the fact 
that Xf = iay\ 

y_(t)W/y'_J(t) = =f 2—^ sin^ cjfJ ± i f 1 — 2 j sin^ cjf j Uy ± z—sin 2a;rcr2. (51) 
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Now, by comparing fl50|) and dST]) . we obtain the following set of equations: 

cos 2na = =p 2 ^ sin^ uf 

< sin 2na cos /3 = ±^1 — 2(^)^ sin^ wr j (52) 

sin 2na sin d = ± - sin 2u:f 

which can be analytically solved in simple cases. The details of the computation when n = 
0,1, 2 are reported in Appendix [B], as well as the values of t and f, which fully characterize 
the control strategy, the corresponding hnal time tf, and the range of existence of these 
extremals in terms of coq and 7 . 

B. Case with Xf = icjx 

In this case cp = |, and, from (155]) . the hnal condition is = bx{0). It is possible to 

prove that a constant control strategy is possible only if wq = 0 , which has been excluded 
at the beginning. Therefore, any control strategy for this family of nonsingular extremals 
must have at least one switching time, and in fact, there must be an even number of them. 
The solutions of ([ 6 |) read 

X,(t) = r+(f)y7f)(y+(t)yi(f-))V+(t), 

XM = r_(i)yt(()(y_(*) 7 (i))V_(i), (53) 

where n = 0,1,2,.... Working as before, and considering the cases Xi{tf) = Xf and 
X 2 {tf) = ^Xf, we hnd the necessary conditions 

K!(f)(u(()ri(i))” = Yl(i)X!Yl(i), 

(y+(*)ri(«))V(() = TV(f')A'a-(*). (54) 

and the total time for the transition is given by tf = 2i + {2n + l)t. By using fl50|) and fl56|l 
we compute 

= — cos {2n + l)a{sin l3I — icos/3ax) — ism{2n + l)aaz, 

Y+{t) = — cos ( 2 ? 7 , + l)a(sin/3/+ z cos/^CTa;) — isin (2n + 1 )q;(T 2 , (55) 
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and the right hand sides of fIMD can be written as 

Yl(t)XfYl{i) = — — sin 2c(;r/+ i — 2 j sin^ ax — sin^ ufa^, 

Y_{t)XfY_{t) = — sin2a;f/+ i — 2 j sin^Ua; + 2i—^sin^wfcr^, (56) 

where h = |. By imposing eqnality of fl55l) and fl56ll . we hnd that Xiitf) = —X 2 {tf) is 
inconsistent, and Xi(tf) = X 2 (tf) = ±^/ is satished nnder the constraints 

' 

cos ( 2 n + l)a sin (3 = sin 2a;f 

< cos(2?7, +1)q!Cos/ 3 = ± — 2 (^)^sin^a;rj (57) 

sin ( 2 n + l)a = ± 2 ^ sin^ cof 

In general, the solntion of this system can be performed nnmerically. Nonetheless, in the 
simple cases n = 0,1 analytical solntions are possible. We refer to Appendix O for the 
details. 

V. SINGULAR EXTREMALS FOR Xf = XswAP 

In this section we consider the case of singnlar extremals, that is, extremal trajectories 

containing at least one singnlar arc. Althongh we limit onr attention to X/ = Xswap, the 

analysis of singnlar extremals can be performed for any hnal target operator. For sake of 

simplicity, this generalization will be presented elsewhere. 

In Section m we have seen that, in the present context, the condition for a singnlar arc is 

bx = by = bz = 0 . On this arc, bx = by = bz = 0 , which implies from flTSll that Cx = Cy = 0 , 

and then Cz = 0, that is, Cz is a non zero constant. It cannot vanish becanse, in the 

Pontryagin Maximnm principle, M 7 ^ 0. Moreover, to satisfy flTSl) . it mnst be Ux = Uy = 0, 

and the only control possibly different from zero is Uz- 

On these extremal trajectories it mnst be L = 0, since this condition holds on any 

embedded singnlar arc, and L is constant everywhere. On singnlar arcs, the Pontryagin 

Hamiltonian gives A = ujoCz, therefore the expression (l2^ is generically valid. We maintain 

2 

onr simplifying assnmption E = ‘^ (which is trne withont loss of generality). 

Now we prove that any singnlar arc must be preceded and followed by nonsingular arcs. 
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Proposition V.l. If Xf is the SWAP operator (or an equivalent operator, according to 
Proposition III. 1 |) . a singular arc cannot be at the beginning, or at the end, of an extremal 
trajectory. 


Proof: Assume that a singular arc is at the beginning of the trajectory. The expression of 
X associated to this class of hnal operators is given in ( 15 ^ . It follows that Cz(tf) = —€2(0). 
Since we have chosen E = ^, it must be Cz{tf) = 1 , and then Cz(tf) = — 1 , which in 
turn implies = 0 . But this is inconsistent with the constancy of the Pontryagin 

Hamiltonian, since its values at t = 0 and t = tf are different. A completely analogous 
argument rules out extremal trajectories ending with singular arcs. 


As we did before, by relying on Proposition III . 31 we can assume without loss of generality 
that b( 0 ) = (63,(0), 0 ), and 63,(0) ^ 0 . The analysis of the dynamics of h{t) on non-singular 
arcs when L = 0 has been detailed in the previous section. We can adapt this treatment 
to the present case, by considering the additional requirement that singular arcs should 
be smoothly connected to nonsingular arcs. In the (63,, 6^) plane, nonsingular arcs follow 
segments through the origin, and singular arcs correspond to a point at rest in the origin. 
Therefore, to connect to a singular arc, the preceding arc must end with 63,(6) = 0 and 
bxit) = 0 . For the hrst singular arc, 63,(6) = 0 gives fITT]) . and 63,(0) = 0 means 

— ojAq sin ijjt± + ojBq cos u}t± = 0 . ( 58 ) 


Joint consideration of these constraints leads to 63,(0) = this is the only initial condition 
such that the system admits singular arcs. Since 63,(0) ^ 1 , we conclude that it is possible to 
have extremals with singular arcs only if 7 ^ wq, which physically corresponds to a control 
power of the same or more magnitude than the interaction. Correspondingly, we have 






Bo = ± — \/T 

70; 


UJ. 


0) 


Cn = 


7^0 
a;2 ^ 


and both the smaller and next to smaller zeros of 63,(6) satisfy 


Ao 

cosa;6± = 

Co 


. , Bo 

sma;6± = —yy- 
Co 


(59) 


(60) 


Without loss of generality we assume that 6+ > 0 , therefore a;(6+ — 6_) = 27 r and the relevant 
times are 6 = 6+, 6 = ^. We can now prove that multiple singular arcs are not possible if 
the extremal trajectory trajectory is optimal. 
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Proposition V.2. A time optimal trajectory may contain at most one singular arc. 

Proof: If a non-singular arc is in between two singular arcs, the corresponding evo¬ 
lution requires a time kr, for some positive integer k. The evolution operators are 
given by Y^{kf) = —I and Yl{kf) = —I, a ccording to flT 6 |l . These contributions are in¬ 
consistent with the requirement of optimality. They amount to no evolution in positive time. 


Propositions IV.II and IV.2I signihcantly constrain the form of any possible singular op¬ 
timal candidates: it must be the succession of evolutions on a nonsingular arc for time t, 
on a singular arc for time t', and hnally on a nonsingular arc, again for time t, because of 
the requirements on the hnal conditions bk{tf) and Ck{tf). In the {bx,by) plane, the state 
starts from the point (^) and move to ( 0 , 0 ) along a straight line, during the hrst nonsin¬ 
gular evolution; then it remains in ( 0 , 0 ) during the singular evolution; hnally, it moves to 
( —^ cos 2(p, ^ sin 2(p) along a straight line, during the second nonsingular evolution. The 
complete evolution on the Lie group reads 

Xi(t) = Z+{i)W+{t')Y+{i), X,{t) = Z_{t)W_{t')Y_{t), (61) 

with 


W+{t) = 

W {t) = 

Z+{t) = 

Z_{t) = 


( 62 ) 


where Uz can be assumed to be constant (see below), and 

/?_!_ VjySy “ 1 “ UJqSz)'} R— UySy (jJqSz^. 

U OJ 

i?_i_ and are functions of the constant controls on the last nonsingular arc, = —7 cos 2 (p 
and Uy = 7 sin 2 (p. 

The evaluation of the two terms in fl6Tl) is cumbersome and the details are reported in 
Appendix [Dl Several extremals are possible, generating any target operator equivalent to 
SWAP. Nonetheless, for the sake of minimizing the transfer time, we can limit our attention 
to the transition to Xf = iay, and compute the corresponding control strategy and transfer 
time tf in this case. 
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VI. OPTIMAL SOLUTIONS FOR Xf = Xswap 


In this section we sum up the previous results and derive the time optimal control to 
reach the SWAP operator for the system. We have found that there is no loss of generality 
considering by = 0 , in which case the only possible target operations are Xf = iay or 
Xf = iaxi ci-nd we have fully characterized the nonsingular and singular extremals, with 
integral of motion L = 0. 

According to Proposition III. 11 to these extremals there correspond extremals for any 
target operator Xf equivalent to the SWAP operator. Therefore, to decide which is the 
optimal solution for some given control strength and drift, it is sufficient to compare the 
hnal times associated to these extremals, and hnd their minimum. We will denote the result 
by topt- Moreover, the cases with Xf = iay and Xf = ia^ correspond to situations in which 
there are an odd, respectively even, number of switching times for the candidate optimal 
control strategies. In the following, we denote the number of switches by s, and use it to 
characterize an extremal trajectory whenever we want to make reference to an arbitrary 
operator equivalent to SWAP. 

First of all we derive lower and upper bounds on topt on nonsingular extremals. 

Theorem 2. On nonsingular extremals, if a control strategy requires s ^ 1 switching times, 
then the optimal time for the transition to Xf = Xswap (or equivalent operator) must satisfy 

(s - 1)- ^ topt ^ 2s-. (64) 

OJ u 

Proof: This result is a direct consequence of Proposition IIV.II For sake of clarity, we 
report the basic facts. If Xf = iay, since tf = 2t + 2nf with n = 0,1,..., it must be 
2n^ ^ tf ^ 2{2n + 1 )^. Since h^ftf) = — 63 ,( 0 ), there is an odd number of switching 
times, that is s = 2n + 1, and fl64l) follows. If Xf = ia^, from tf = 2t P (2n + l)t we hnd 
( 2 n + 1 )^ ^ t/ ^ 4 (? 7 , + 1 )A Since bx{tf) = 63 ,( 0 ), in this case there is an even number of 
switching times, s = 2(n + 1), and fIMll is still valid. 

We have found that there are not extremal trajectories without switching times, s = 0. 
We have analytically derived the extremals (and corresponding control strategies) when 
s ^ 5. We limit our analysis to these cases, and, in view of Theorem |2l we can hnd the 
optimal solution if tf ^ which is the lower bound for topt if s ^ 6. In fact, under this 
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tf topt 



FIG. 1: Left plot: comparison of the transition times tf for the extremal trajectories derived in 
this work, as a function of 7 , when ujq = \. The blue and black lines are associated with the cases 
of even, respectively odd number of switching times in nonsingular extremal strategies, and the 
red line corresponds to the singular extremal. The dashed lines represents some upper and lower 
bounds on tapt as derived in Proposition [2j More precisely, they give for k = 2,3,4, 5 and 6 . 
Right plot: resulting optimal time in the range of values of 7 where analytical expressions for the 
extremals are possible. 


condition, the extremal trajectory that we have not analyzed cannot outperform those which 
we have taken into account. 

It is not possible to give the form of the time optimal control and minimum time in a 
compact way since the optimal solution depends in an intricate way on the relative values of 
Uo and 7 . In different situations the optimal solution is given by different types of extremals, 
singular or nonsingular. This is illustrated in Fig. [H where the choice Wq = 1 has been made, 
and the extremals which are clearly sub-optimal have been neglected. We summarize our 
hndings in the following theorem, where, on account of Proposition III. 11 we provide the 
result for any operator equivalent to SWAP. 


Theorem 3. Consider the system Given a target operator of the form 






(65) 


for some 9 ? G M, the optimal control strategy steering Xi and X 2 to can be nonsingular 
or singular, depending on the relative values of uq and 7 . 
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Nonsingular strategies are hang-hang, with s switehing times. They are given hy 

tta; = ± 7 cos 99 , tiy = =F 7 sin 99 , = 0 , ( 66 ) 

and both Ux, Uy ehange sign at tk = t kt, with k = The number s depends on 

the relative values of uq and 7 . The optimal transition time is given by topt = 2t + (s — l)f 
for a given s. Analytieal expressions oft and t are provided in AppendicesWl andllA 

The optimal singular extremal eonsists of the alternation of a nonsingular arc for time t, a 
singular are for time t', and finally a nonsingular arc for time t. The optimal transition time 
is given by topt = + t', and analytical expressions oft and t' can he found in Appendixl^ 
The associated control strategy is given hy TUB) on nonsingular arcs, and Ux = Uy = Uz = 0 
on the singular arc. 


VII. DISCUSSION AND CONCLUSIONS 

The problem of controlling in minimum time a quantum bit interacting with an additional 
two level system can be mapped to an equivalent problem of simultaneous control of two 
quantum bits under the same control held. In this paper, we have solved this problem for 
the class of SWAP operators, although some of our results apply to the case of an arbitrary 
hnal target operations as well. The application of the necessary conditions of optimality, 
given by the Pontryagin Maximum Principle, leads to the consideration of two type of 
candidates optimal controls, singular and nonsingular. Our analysis shows that singular 
extremals can be optimal only in a specihc range of control strength, and they must have 
a particular structure, with a singular arc between two nonsingular arcs. The dynamical 
analysis benehts from known results from the central force problem. We have found that 
the optimal control strategy depends in a complicated way on the relative strength of the 
control and the interaction. These results can be used to synthesize the SWAP operation in 
minimum time in quantum computation implementations in cases where the target system 
is interacting with an analogous system, and-or to synthesize the SWAP operation on two 
systems simultaneously. 

Because of the homomorphism between SU{2) and SO (3), we can adapt our results to 
the derivation of optimal control strategies for rotating a rigid body in in minimum 
time. The class of SWAP operators in SU{2) corresponds to the class of rotations of angle 
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TT around axes orthogonal to the axis of the drift dynamics. This is another scenario where 
our analysis applies. 

The specihc form of the interaction (Ising) is crucial for the derivation of our results, 
and for the equivalence of the two problems mentioned above. An interaction term of the 
form cTi ® (T 2 , with arbitrary a\ and cr 2 , can be considered under the perspective of driving 
only one two-level system, the second one representing an undesired disturbance. But in 
this case this problem is not equivalent to that of simultaneously driving a pair of two- 
level systems with opposite drifts. Finally, more complicated interaction terms (e.g., the 
Heisenberg interaction) are not compatible with our procedure. 

In our analysis, the drift parameter wq and the control strength 7 can be arbitrarily varied, 
and we have analytically derived the optimal control strategy and the corresponding time 
when — exceeds a threshold, which is approximatively 0.325. For smaller control strength. 


numerical investigations are needed. 

We have found that the optimal control strategy for the SWAP operator (or equivalent) 
always requires = 0, both for nonsingular or singular arcs. Therefore, we can conclude 
that, for this class of operators, the scenario with only two independent controls Ux and Uy is 
completely equivalent to that considered in this work. The same result has been found when 


B. 


and 


investigating the analogous optimal control problem with only one two-level system 
it is a consequence of the symmetry of the problem. 

Our results add to the growing literature on optimal control of quantum systems. The 
problem of time optimal simultaneous control of two quantum bits in minimum time does 
not seem to have been considered earlier. One notable exception is 28|], where however the 
problem was set up for the state and not for the evolution operator as here. 


[1] F. Albertini and D. D’Dalessandro, J. Math. Phys. 56, 012106 (2015) 

[2] R. Romano, Phys. Rev. A 90, 062302 (2014) 

[3] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cam¬ 
bridge University Press, Cambridge, U.K., New York (2000) 

[4] M. H. Levitt, Spin dynamics: basics of nuclear magnetic resonance, John Wiley and sons. 
New York-London-Sydney (2008) 


24 




[5] R. Wu, C. Li and Y. Wang, Phys. Lett. A 295, 20 (2002) 

[6] U. Boscain and P. Mason, J. Math. Phys. 47, 062101 (2006) 

[7] M. Wenin and W. Potz, Phys. Rev. A 74, 022319 (2006) 

[8] A. Carlini, A. Hosoya, T. Koike and Y. Okudaira, Phys. Rev. A 75, 042308 (2007) 

[9] E. Kirillova, T. Hoch and K. Spindler, WSEAS Trans. Math. 7, 687 (2008) 

[10] A. Garon, S. J. Glaser and D. Sugny, Phys. Rev. A 88 043422 (2013) 

[11] G. G. Hegerfeldt, Phys. Rev. Lett. Ill, 260501 (2013) 

[12] G.D. Aiello, M. Allegra, B. Hemmerling, X. Wang and Paola Gappellaro, arXiv: 1410.4975 
(2014) 

[13] B. Russell and S. Stepney, J. of Phys. A 48, 115303 (2015) 

[14] D.A. Lidar, Adv. Ghem. Phys. 154, 295 (2014) 

[15] G. de Lange, T. van der Sar, M.S. Blok, Z.H. Wang, V.V. Dobrovitski and R. Hanson, Scientific 
Reports 2, 382 (2012) 

[16] T. H. Taminiau, J.J.T. Wagenaar, T. van der Sar, E. Jelezko, V.V. Dobrovitski and R. Hanson, 
Phys. Rev. Lett. 109, 137602 (2012) 

[17] G. Arenz, G. Gualdi and D. Burgarth, New J. Phys. 16, 065023 (2014) 

[18] L. Pontryagin et V. G. Boltyanskii, Mathematical theory of optimal processes, Mir, Moscou 
(1974) 

[19] W. Eleming and R. Rishel, Deterministic and Stochastic Optimal Control, Applications of 
Mathematics, Springer-Verlag, New York, 1975 

[20] U. Boscain and B. Piccoli, Optimal Syntheses for Control Systems on 2-D Manifolds, Springer 
SMAI 43 (2004) 

[21] R. Wu, R. Long, J. Dominy, T.S. Ho and H. Rabitz, Phys. Rev. A 86, 013405 (2012) 

[22] M. Lapert, Y. Zhang, M. Braun, S. J. Glaser and D. Sugny, Phys. Rev. Lett. 104, 083001 

( 2010 ) 

[23] R. Broucke, Astrophisics and Space Science 72, 33 (1980) 

[24] D. D’Alessandro, Introduction to Quantum Control and Dynamics, CRC Press, Boca Raton 
EL (2007) 

[25] N. Khaneja, R. Brockett and S.J. Glaser, Phys. Rev. A 63, 032308 (2001) 

[26] D. DAlessandro and M. Dahleh, IEEE Trans. A. C. 46, 866 (2001) 

[27] R. Huneault, Time Optimal Control of Closed Quantum Systems, Master thesis, Department 


25 


of Mathematics, University of Waterloo, Ontario, Canada, 2009 
[28] E. Assemat, M. Lapert, Y. Zhang, M. Braun, S.J. Glaser and D. Sugny, Phys. Rev. A 82, 
013415 (2010) 


[29] 


[30] 

[31] 


[32] 

[33] 


Although we assume that three independent control actions can affect the system’s dynamics, 
we will find that our results are valid also in the case with only two independent control actions 
entering the dynamics, when we limit our attention to Xf = XswAP- 

The case of a final operation Xf gU(2) can be dealt with by adapting the results in SU{2) 
For the unbounded minimum time control problem, one can find the optimal control in the 
sense of 25|. The form of the optimal control can be significantly different for different types 


of bounds on the control such as bounds on the Li or Lqo norm. 

The multiplicative constants are introduced for further reference. 

Alternatively, we could study the costate dynamics by using the polar coordinates b and 9. 
While this approach makes more transparent the symmetry of the problem (which is, basically, 
independence from 9), we prefer to present our analysis in terms of bx and by, because they 
have a smooth evolution. This is not the case for b (whose derivative changes sign when 6 = 0) 
and 9 (which flips between two fixed values). 


Appendix A: Numerical analysis supporting Conjecture 1 


In this appendix we provide numerical evidence that L = 0 is a necessary condition for the 
optimality of extremal trajectories. Without loss of generality, we write the initial condition 
as (6a:(0),0), which, with the choice E = leads to 

bx{0) = cosd\/l - 6x(0)2, by{0) = sin^?^l - bx{0y (Al) 


for some "d satisfying 0 ^ d ^ tt. We solve the costate dynamics (j26jl . and from them we 
derive the corresponding extremal controls (El])- Using the extremal controls, we numerically 
solve Eq.s (|6|) and derive the corresponding extremal trajectories in the Lie group, Xi = 
Xi(t,'9,bx(0)) and X 2 = X2{t,'i},bx{0)), in the time interval 0 ^ ^ T. We then consider 

the following real functions: 


F. 


+ 


F_ 


F+{t,^,bxm 

F.{t,^,bxm 


1 

4 

1 

4 


(Tr(Xi - X2)ax) 
(Tr(Ai + X2)ax) 


2 

2 


+ (Tr(Xi-A2)a,)'], 
+ (Tr(Xi+X2)a,)'], 


(A2) 
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satisfying 0 ^ ^ 1. Moreover, the function F+ reaches its maximum if and only if 

Xi = X2 matches some operator equivalent to SWAP (according to Proposition lll.ljl . and 
similarly for F_, with Xi = —X 2 . Therefore the optimal time tf we want to determine is the 
smallest time such that F+ = 1 or F_ = 1, and it can be found by gradually increasing the 
range of numerical integration [0,T]. For any choice of ojq and 7 that we have considered, 
we have always found that the smallest time is associated with'd = 0 or tt, or 6a;(0) = 0 or 
1, which means L = 0. For sake of completeness, some contour plots of F± are reported in 
Fig.s[6]and[71 We do not numerically compute the minimum time tf, because once we know 
that L = 0 we can obtain an analytical result. 

This numerical analysis is fully consistent with the analytical expressions of nonsingu¬ 
lar extremals presented in Section IIVI With suitable choices of T, that is, by sufficiently 
increasing it, it is possible to investigate the structure of sub-optimal extremal arcs, and, 
again, we hnd consistency with our analytical results (see Fig. [7]). On one side our numerical 
computations motivate the analytical evaluation of extremal strategies; on the other, they 
provide an independent check of these results, since they are fully consistent with them. 

Appendix B: Analysis of nonsingular extremals for Xj = iay 

System fl52|) admits analytical solutions for several values of n. We recall that a necessary 
condition for all these extremals is X 2 {tf) = —Xi{tf). 

The case n = 0 is the only possible extremal with only one switching time for u^- We 
hnd that system 0521) is inconsistent, with the only exception of the case when 7 = wq- In 
this special case, it must be coswf = 0 and Xi{tf) = —Xf, and we hnd tf = 2i = \/2^ and 
6a;(0) is unconstrained. 

If n = 1, the extremal control strategies have 3 switching times. We hnd that 



(Bl) 


where the signs of the two terms follows from Proposition lIV.il which implies that sin ojt ^ 0, 


cos UT ^ 0, and sinwr ^ 0. The upper signs correspond to Xi(tf) = Xf, the lower signs to 
Xi{tf) = —Xf. Eq.s OBip are consistent with flTT|) and fH3|) when 
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FIG. 2: Plots of ^^(O) and as functions of the control strength 7 , with ojq = 1, Xf = iay and 
n = 1. Dashed lines are associated with extremals which are certainly sub-optimal. 


These solutions exist only when \/2 — 1 ^ ^ ^ 1, if we consider the upper sign, and when 
1 ^ ^ ^ \/2 + 1 if we consider the lower sign (otherwise the two functions in fIBip are 
not dehned, or bx{0) < 0). The total time for the transition is given by tf = 2{t + t). For 
Xi{tf) = Xf we hnd 
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Similarly, for the transition Xi{tf) = —Xf we obtain 
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t/wq + 7 ^ V 2wo V 7 '^7 \ ^0 

In fIBSp and flB4p . we have nsed the fact that the sign of coswf is determined by the third 
eqnation of ([52]): it is positive in the hrst case, negative in the second case. Plots of 63 ; (0) 
and tf as fnnctions of 7 are shown in Fig. | 2 j 

If n = 2, the extremal control strategies have 5 switching times. By expressing cos4a: 
and sin 4a in terms of sin a and cos a, we rewrite system (j52ll as 


8 sin^ a (l 


sin^ a 
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sm a] 


4 sin a (1 — 2 sin^ a) cos wr = =F^ sin 2wr 
From the hrst and second eqnation we derive the biqnadratic eqnation in sin a. 
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FIG. 3: Plots of ^^(O) and as functions of the control strength 7 , with ojq = 1, Xf = iay and 
n = 2. Dashed lines are associated with extremals which are certainly sub-optimal. 


which admits four solutions. If Xi{tf) = 

sin a = 


X 

\\ 


f, we obtain 
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Conversely, if Xi{tf) = —Xf, we hnd 

From Eq. fHHj) we can find sin uif, and moreover, from the first equation in fIBSp , by choosing 
the suitable signs, and rewriting the left hand side as cos da for simplicity, we can evaluate 
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■ cos da. 


These conditions are consistent with fjdTl) and fjd3ll as long as 


ftx(O) = 


,2 0^2 , 


7^ cos^ a 
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(BIO) 


Now, by requiring 0 ^ sin a ^ 1 and 0 ^ ^x(O) ^ 1, it is possible to prove that all the 
solutions of (IBOp are consistent, and to determine the associated ranges of which are 
not reported here because they are not needed for further developments. The hnal times 
associated to these four extremals are given hj tf = 2i + df, that is 
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Plots of 63 ;( 0 ) and f/ for these four extremal trajectories are shown in Fig.|3l Our development 
shows that the only extremal which can contribute to optimal solutions corresponds to the 
solution with positive sign in Eq. fIBSp . which is associated with Xi{tf) = —Xj. 


Appendix C: Analysis of nonsingnlar extremals for Xf = iax 


We compute the analytical solution of system flFTp in the simplest cases. We recall that 
in this case it must be Xi{tf) = X 2 {tf) = Xf. 

If n = 0, the extremal control strategies have 2 switching times. By considering the 
dehnitions fHHl) . the system is rewritten as 

cos CUT = =F-sin2a;f 

' U) 

< ^sinwr = ± (l-2(^)'sin2a;f) . (Cl) 

^^sinwf = ±2^sin2a;f 

Consider the case of upper signs in (IClll . that is, hnal operator Xi(tf) = Xf. Since ojq 7 ^ 0, 
we hnd sinwr = and sinwr = which are well dehned only when ^ ^ 1, or 
^ Moreover, from sinwr = ± sin cur and the hrst equation in flClD . we hnd that 
cos wf = =F cos Ljf. We conclude that cos ui = cos ut and sin u)i = — sin ut. These are 
constraints on the switching times, which must be consistent with the expressions fj4Tl) and 
fHSp . By using the explicit expressions on Aq, Bq and Cq, we hnd that the initial condition 
of the costate must satisfy 


bx{0) 


^ / (n;g-27^)^-7~ 
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Therefore, this extremal is well dehned when 7 ^ and the total time for the transition 

is given hj tf = 2t +1, or more explicitly 

2 f 

tf = — , TT + arctan 

The case with hnal operator Xi{tf) = —Xf is obtained by replacing r with r + tt in the 
previous formulas. However, the result sinwf = — ^ < 0 is inconsistent with the conditions 
expressed in Proposition IIV.II therefore there are no extremal trajectories in this case. See 
Fig. mfor plots of 63 ,( 0 ) and t/ in this case. 


^o + T 

372 - 


(C3) 


30 















bAQ) tf 




FIG. 4: Plots of ^^(O) and tf as functions of the control strength 7, with wq = 1, Xf = iax and 
n = 0. Dashed lines are associated with extremals which are certainly sub-optimal. 


If n = 1, the extremal control strategies have 4 switching times, and the system (l57|l 
becomes 

(4 sin^ a — 1 ) cos wf = sin 2uf 

< ^ (1 — 4sin^ a) sin a = ± — 2 sin^ turj . (C4) 

sin a (3 — 4 sin^ a) = ±2^ sin^ uf 

From the second and third equation of this system we get a depressed cubic equation in 
sin a, 

sin^ a-f-psina-|-g = 0, (C5) 


where 


P= 


q = ± 


Uo 


(C6) 


2’ ^ 87 

This equation can be solved by using the Cardano’s method. By dehning u and v such that 
M -f v = sin a, we hnd = —q and 3uv = —p, we hnd that both and satisfy the 

equation 


z +qz- 


3/ 


= 0 . 


(C7) 


Then we solve this quadratic equation and get and v^, extract the cubic roots, and 
impose the further requirement 3uv = —p. We are interested in real roots of the original 
cubic equation, and their number depend on the sign of the discriminant A = —27g^ — 4p^. 
If A < 0, there is only one real root; if A = 0 there are three real roots with degeneration, 
and hnally if A > 0 there are three distinct real roots. By considering the values of p and q 
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in (1C6|1 . we conclude that, if ^ < fyf, the only real root is given by 


sin a = 



(C8) 


where we consider the real cubic root. Since it must be sin a ^ 0, we require g ^ 0, and then 
the only possibility is Xiitf) = —Xf, corresponding to the lower sign in flC4p and flC6p . If 




- 1 / |, all the roots are real, and their general expressions are 


{sma)k = 2^J-^cos (^6 (C9) 

with A; = 0,1, 2, and 

, 1 1 
0 = — arctan — 

3 3q 

In the given range of parameters, sin Set is positive for A: = 0 and negative for k = 1,2. 
Therefore, in the hrst case we have to consider the upper signs in flC4p and flCbp . leading to 
Xi{tf) = Xf, and the lower signs otherwise, with Xi{tf) = —Xf. 

The requirements 0 ^ sin a ^ 1 and 0 ^ sin ojf ^ 1 further restrict the range of admissible 
values for Their explicit expression are cumbersome and unnecessary; we only mention 
that extremal trajectories associated to the case A: = 0 in (IC9D are impossible because of 
the constraint sinwf ^ 1. Therefore, Xi{tf) = —Xf is the only possible hnal operation. By 
solving the last equation in (IC4p we hnd 



sin cur = 




2uoj 


sin 3a;, 


(Cll) 


which is well dehned because sin 3a ^ 0. The conditions fICSp and fICllD are consistent with 
flTT]) and fl43p as long as flBlOp is satished. The total time for the transition on this extremal 
trajectory is given hy tf = 2i + 3i, that is 
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37 r + 2 arcsin 




2uoy 
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3 arcsin ( — sin a 
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Fig. E] shows plots of &a;(0) and tf for this class of extremals. According to the analysis of 
Section ED the only extremal which can be optimal is the one characterized by 
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FIG. 5: Plots of ^^(O) and tf as functions of the control strength 7 , with wq = 1, Xf = iax and 
n = 1. Dashed lines are associated with extremals which are certainly sub-optimal. 

Appendix D: Analysis of singular extremals 


Considering the two terms in fl 6 T|) . and assnming that Bo ^ 0, we can be recast the hnal 
resnlt in a rather compact form as 

Xi{t) = — ^sinr^sin -f r;) — cos(pcos -|-(p)j/-|- 

- ^ ^(1 - cos 2 (f) cos (^+ + 17) + sin 2 (f sin (,^+ + r])^ax + 

- ^ ((1 + cos 2ip) sin (,^+ + v) + sin 2ip cos (,^+ -b 17 ) j ay + 

- z ^ sin 77 cos (^+ + 77) — sin f cos (^+ + f)^ az 


and 


X 2 it) = -|-^sin 77 sin — 77 ) + cos</ 9 cos (^_-b -b 

— ^ (^(1 “ cos 2(f) cos ((^_ — 77 ) + sin 2f sin (^_ — 77 ) j 0 - 3 ; -b 

— ^(1 -b cos 2f) sin — 77 ) -b sin 2f cos — 77 )j ay + 
-b 7 ^ sin 77 cos — 77 ) -b sin (f cos (^_ -f 72 ) j a^ 

where we have dehned the two fnnctions 
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Since Uz enters the problem only through the integral in fIDSp . without loss of generality 
we can assume it is constant on the singular arc. Therefore = ^(uz ± uJo)t'. The free 
parameters in flDip and flD2p are (p and (or, which is the same, (p, Uz and t'). Target 
operators equivalent to SWAP have the form 

Xf = = i{siinpax + coscpcTy), (D5) 

and, by requiring Xi = ±Xf we hnd the system 

sin rj sin (,^ 4 . + //) — cos (p cos (^+ + <^ 9 ) = 0 

(1 — cos 2 (p) cos ((^+ + v) + sin 2(p sin ((^+ + h) = T 2 sin (p 
< (D 6 ) 

(1 + cos 2ip) sin (^_|_ + rj) + sin 2ip cos (^4 + 77 ) = =f2 cos ip 

^ sin r] cos ((^4 + r]) — sin ip cos (^4 + (p) = 0 

From the hrst and fourth equations we hnd that cos (^4 + p>) = ± sin 77 , and consequently 
sin (,^4 + 77 ) = ±cos(p and cos (^4 + 77 ) = ±sin(p. We observe that the second and third 
equations in (1D6H reduce to identities. Finally, the solution of (ID 6 P is expressed as 

cos ^4 = ± sin ((79 + 77 ), sin ^4 = ± cos ((p + 77 ). (D7) 

For solving X 2 = iW/, we follow the same steps. We hnd that 

cos^_ = =F sin ((p — 77 ), sin^_ = =F cos ((p — 77 ). (D 8 ) 

Since the signs in (lD7p and (1D811 are opposite, the only consistent solutions must satisfy 
Xi = —X 2 = Xf, with Aj as in fIDSp for some (p. By combining Eq.s flD7p and flD 8 p . we hnd 
that Uzt' = 2tt — 2pi-\-Akir or Uzt' = An — 2ip + Akn for some integer k, and WqF = n — 2ri + Aln 
or uqI' = 3n — 21] + Ain for some integer / ^ 0. We are interested in the solution with 
minimum time t', with the additional constraint that \uz\ ^ 7 . This minimum is obtained 
when (p = 0, A: = —1 and I = 0, that is, ujot' = n — 2r] and Uz = 0. This means that, on 
singular arcs on extremal trajectories, the evolution is given by the drift term. 

Summing up, the total time for the transition X ^ Xf = i(jy is given hy tf = 2t + t', 
where t is determined by fl60l) . The explicit expression in terms of Uq and 7 is 

tf = — (n — arccos f^ ) H- (n — 2 arcsin —^ . (D9) 

' 77+7V \i) ) ^<>\ '1) 
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When i?o > 0 , the analysis is completely analogous. The dynamics in the Lie group is 
described by the following equations: 

= T^sinr^sin (^+- p) + cos(pcos (^++ J + 

+ ^ (^(1 - cos 299) cos (^+ -r]) + sin 2 (p sin (^+ - ??) j + 

+ ^ (^(1 + cos 2ip) sin ((^+ - Tj) + sin 2ip cos (^+ - p) j ay + 

+ f ^sinr7cos ((C_,_ — 77) + sin(^cos (.^++ 99)ja'2 (DIO) 


and 


^•2if) = — ^sinpsin ((^_ + 77) — cos(pcos ((C_ + 7?)^/+ 

+ ^ (^(1 - cos 2 ip) cos (^_ + 77) + sin 2 (p sin + 77) j a^ + 

+ ^ ((1 + cos 2 ip) sin + 77) + sin 2 ip cos + 77) j ay + 

— 7^ sin 77 cos (^_ + 77) — sin 7?cos (^_ + (p^az- 


(Dll) 


As before, the only consistent requirement is Xi = —X 2 = Xf, and the solution is given by 


cos ^_ = ± sin (79 + 77) , sin (C_ = ± cos (79 + 77). 


(D12) 


For solving X2 = ^Xf, we follow the same steps. We hnd that 

cos^+= =F sin (79 — 77), sin^+= =1= cos (79 — 77). (D 13 ) 

By combining these equations, we find that Uzt' = 2 tt — 2ip Akn or Uzt' = dyr — 279 + Akn 
for some integer /c, and ojQt' = vr + 277 + Ain or uj^t' = Stt + 277 + Ain for some integer 1^0. 
These extremals are certainly not optimal, as both t and t' are larger than the corresponding 
quantities derived in the case Bq ^ 0. 


35 








(9 


F. 



IS 




FIG. 6 : Contour plots of the functions F+ and F- for selected values of ujo, 7 and T (in the inset). 
The dashed lines represent points with a fixed value for L, decreasing towards the border of the 
plot, where L = 0. The red spots shows the values of 5a;(0) and such that F± = 1 with 0 ^ t ^ T. 
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FIG. 7: Contour plots of the functions F+ and F_ for selected values of ujq, 7 and T (in the inset). 
The parameters in the two plots are the same, but T which is increased. From the behavior of FT, 
it is apparent the existence of a sub-optimal extremal. 
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